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Quantum Chromo Dynamics (QCD), the theory of strong interactions, predicts a transition of the usual matter to a new 
phase of matter, caUed Quark-Gluon Plasma (QGP), at sufficiently high temperatures. The non-perturbative technique of 
defining a theory on a space-time lattice has been used to obtain this and other predictions about the nature of QGP. Heavy 
ion collisions at the Large Hadron Collider in CERN can potentially test these predictions and thereby test our theoretical 
understanding of confinement. This brief review aims at providing a glimpse of both these aspects of QGP. 



1. Introduction 

There are two very commonly quoted motivations for 
the upcoming Large Hadron Collider (LHC) at CERN 
in Geneva, the center of attraction for the articles in 
this volume. Perhaps the primary one is that LHC wiU 
provide us a key to understand the origin of the visible 
mass of our Universe. This alludes to that fact that 
our standard model(SM) of particle interactions has to 
start with matter in the form of massless quarks and 
leptons. The famous Higgs mechanism [ 1 of sponta- 
neous breaking of gauge symmetries provides masses to 
them, and the carriers of the weak force, namely ,Z. 
LHC is widely expected to discover the Higgs boson 
which is tied with this mechanism. The other moti- 
vation rests on the fact that the standard model has 
been well understood due to the many impressive pre- 
cision tests carried out in many experiments, including 
those at the Large Electron Positron (LEP) at CERN 
and the Tevatron at the Fcrmilab in the USA. However, 
new physics beyond the standard model (BSM) has to 
exist [[5] since SM contains many, at least 19, arbitrary 
parameters and thus cannot be the final theory. Indeed, 
it is even hoped that LHC may provide us not only a 
glimpse of the BSM physics, but it will hopefully also 
explain the origin of the mass of the dark matter in the 
Universe. 

While these motivations are largely correct, there 
are certain oversimplifications in them, leading to a 
few misconceptions, especially in the popular media. 
First of all, even if the expected Higgs particle is actu- 
ally discovered, the origin of the mass of up/down(u/(i) 
quarks can be claimed to be understood only after it 
is also established that the Higgs particle couples to 
them with a strength of 10~^, not an easily achiev- 
able goal at LHC. Indeed, one may as well need an 
electron-positron collider to establish this in the post- 



LHC era. Moreover, the protons and neutrons, which 
make up most of the visible mass in our Universe, have 
each a much larger, almost a factor of 100 larger, mass 
than the sum of the masses of their constituent u/d 
quarks. Therefore, the understanding of the visible 
mass of the Universe will emerge from the efforts to 
figure out why protons/neutrons have such large bind- 
ing energies. Starting from molecules to atoms and 
nuclei, we are accustomed to the idea that the inter- 
actions which bind the respective constituents give rise 
to binding energies much smaller, less than even a per 
cent. This has given rise to the very successful idea 
of treating these interactions perturbatively as an ex- 
pansion in the strength of the interaction. As we shall 
see below, one needs new suitable techniques to investi- 
gate these large binding energies, in Quantum Chromo 
Dynamics (QCD), the theory of interactions of quarks 
with gluons, the carriers of the strong force. 

As may be seen from the articles in this volume it- 
self, QCD is an integral part of our standard model of 
particle and their interactions. From various experi- 
ments in the past, it is well known that quarks carry 
both flavour quantum numbers such as, electric charge 
or strangeness, as well as colour: they transform as a 
triplet under the colour SU{3) group. As in the case 
of electric charge, the colour charge is also mediated 
by massless vector particles, gluons. Structurally, the 
theory of quark-gluon interactions, QCD, looks very 
similar to that of electron-photon interactions, QED. 
A key difference though is that there are eight gluons 
which themselves carry colour charge, transforming as 
an octet under SU (3)-colour group. Consequently, glu- 
ons can interact amongst themselves. Furthermore, the 
QED coupling is rather small at the scales we probe, 
being 1/137, whereas the smallest measured QCD cou- 
pling, as, is about 0.12. In fact, more often, one has to 
deal with — 0.3 or so and it is 1 in the bound 
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states like proton or neutron. QCD exhibits a much 
richer structure and a variety of phenomena as a result 
of this large a^. Quark confinement and dynamical 
chiral symmetry breaking can be named as typical ex- 
amples. A lack of observation of quarks in experiments 
led to the hypothesis that quarks are permanently con- 
fined in the hadrons, i.e, protons or pions whereas the 
lightness of pions compared to protons is expected to 
be understood as the phenomenon of dynamical break- 
ing of the chiral symmetry by the vacuum. QCD as the 
theory of strong interactions has to explain these phe- 
nomena. Since, QCD is too complex, simple models 
based on underlying symmetries are often employed to 
account for its non-perturbative aspects. Indeed, most, 
if not all, of the "precision tests" are either performed 
experimentally only at small coupling, a^., correspond- 
ing to rather rare events, or employ the simple QCD- 
based models. The latter are in many cases possible 
weak links in the precision tests of the standard model 
: physics beyond standard model may even show up 
in non-perturbative QCD beyond these models. We 
need to look for it and rule out such a mundane pos- 
sibility for BSM-physics in order to be sure that other 
exotic possibilities are indeed worth looking for. Thus, 
non-perturbative techniques are needed for real preci- 
sion tests of QCD. As a glaring example, let me men- 
tion that the easiest precise measurement at LHC will 
perhaps be the total proton-proton cross section at 14 
TeV. The current best theoretical prediction for it is 
[[3] cr*°* = 125± 25 mb ! As explained in [[g, one 
uses the so-called Regge Models to arrive at it, and one 
such model can even explain the currently observed Q^- 
variation of the structure function of proton, F2, as well. 
Recall that a key cornerstone for establishing QCD as 
the theory of strong interaction is this variation. 

While obtaining a reliable prediction for the above 
cross section from QCD still seems far away, a non- 
perturbative technique does exist today to obtain other 
quantities, such as the decay constants or the weak 
matrix elements, from QCD using first principles, and 
these could still provide non-perturbative precision 
tests of the standard model. QCD defined on a space- 
time lattice is such a tool. Not only does it explain 
many of the above mentioned phenomena but it pro- 
vides quantitative estimates of many physical observ- 
ables. Furthermore, the same techniques of lattice 
QCD lead to spectacular predictions for the behaviour 
of matter under extreme conditions. Thus, lattice QCD 
predicts the existence of a new phase, called Quark- 
Gluon Plasma (QGP) at sufficiently high temperature, 
and a phase transition of the strongly interacting mat- 
ter of protons, neutrons and pions to the new phase 



QGP at high enough temperature. The dynamically 
broken chiral symmetry of QCD at low temperatures 
in our world is expected to be restored in the QGP 
phase, 'melting' away the constituent mass of the light 
quarks acquired due to interactions. 

Our Universe ought to have existed in such a phase 
a few microseconds after the Big Bang, and about 20 
microseconds later the phase transition to the normal 
hadrons like protons, neutrons and pions ought to have 
taken place in it. Whether there arc any imprints of 
this phase transition on the astronomical objects ob- 
served today depends on the nature of the phase transi- 
tion. There have been speculations of stars with strange 
matter, consisting of neutral baryons made from an up, 
down and a strange quark each. Similarly attempts 
have been made to study the influence of such a phase 
transition on the Big Bang Nucleosynthesis. More ex- 
citingly, the Large Hadron Collider (LHC) itself will 
provide us with an opportunity to create these Early 
Universe-like conditions of high energy densities, or 
equivalently high temperature, in the laboratory in its 
proposed heavy ion collisions of Lead on Lead at 5.5 
TeV colliding energy. Heavy ion collisions at relativis- 
tically high energy have had an illustrious past, and 
even more impactful present. Early such experiments 
were made at the SPS collider in CERN, Geneva at a 
colliding energy of 17 GeV per nucleon in the center 
of mass (cm) frame. The relativistic heavy ion collider 
(RHIC) has been operative in BNL, New York, since 
a last few years and has produced heavy ion collision 
data for a variety of ions. Deuterium (D), Copper(Cu), 
and Gold (Au), at a spectrum of energies, 62-200 GeV 
per nucleon in the cm frame. Experiments at LHC will 
thus see a further jump in the colliding energy by a 
factor of about 30. It is hoped that this will offer us 
cleanest environment yet for investigating the physics 
of quark-gluon plasma. 

In this short review, we shall attempt to provide a 
glimpse of how lattice QCD leads to QGP and predicts 
many of its properties as well as those of the corre- 
sponding phase transition and how the heavy ion colli- 
sion experiments amazingly provide us an opportunity 
to produce QGP in a laboratory, including the expec- 
tations of what we may observe at LHC. 

2. QGP from Lattice QCD 

In order to understand and appreciate the funda- 
mental importance of attempts to discover QGP at 
the LHC, let us first review the basics of lattice QCD 
and why it facilitates a truly reliable treatment of non- 
perturbative physics. In the process, we shall also see 
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why essentially the same tested technique for obtaining, 
say, the hadron masses, conies into play for predicting 
new phases or phase transitions. 

2.1. Basic Lattice QCD 
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Figure 1. Quark and gluon fields on a space-time lat- 
tice. 

Lattice field theory is defined by discretizing the 
space-time. The (inverse of the) lattice spacing a acts 
as the ultra-violet cut-off needed to tame the diver- 
gences in a quantum field theory. One places the anti- 
commuting quark fields ijjix)^ and ip{x) on lattice sites 
whereas the gluon fields reside on the links, as shown 
in Figure [1] A directed link from site x in the posi- 
tive direction /t is associated with the gluon field U^, 
while the link to the site x — p, in the opposite direction 
is U^_p^. A gauge transformation Vx € SU{2>) rotates 
the quark field in the colour space : i^^x) = Vxi^ix). 
Demanding that the gluon field at the link x in the di- 
rection /i, Uy,{x), change to U'^ix) = VxU^,{x)V,^_^f^, en- 
sures that the (discrete) kinetic energy term of quarks 
remains invariant under such a gauge transformation. 
Constructing gauge actions from closed Wilson loops of 
the links, like e.g., the smallest square loop, called pla- 
quette and displayed in Figure [Tl ensures their gauge 
invariance. 

It turns out that a straightforward discretization of 
the derivative, given by [a ■ dfj_ip{x) — ^p{x + afi) — 'ip{x — 
a/t)], can be made gauge invariant as shown in the Fig- 
ure [l] where the links end on respective quark fields ip 
at the sites. Thus a sum over all independent terms of 
both types shown in Figure [1] yields the QCD action on 
the lattice. However, it leads to the so-called Fermion 
Doubling problem : each lattice fermion corresponds 



to 2"^ ~ 16 flavours in the continuum limit of a — > 0. 
Various lattice Fermion actions, referred to as the Stag- 
gered, Wilson, Domain Wall or Overlap Fermions, have 
been proposed to alleviate this problem. In view of 
their simplicity and an exact chiral symmetry even on 
the lattice, the staggered Fermions have dominated the 
field of interest for this article, namely lattice QCD at 
finite temperature and density. Briefly, these are sin- 
gle component Grassmann variables on each site, with 
the 7- matrices replaced by suitably defined sign factors. 
They have at/(l)xC/(l) chiral symmetry and 4 flavours 
in continuum limit. An oft-discussed problem of the 
staggered Fermions, though, is that two or three light 
flavours are not simple to define, and the currently used 
methods may miss out on important physics aspects re- 
lated to anomalies. It is often argued that for the bulk 
thermodynamic properties these issues are likely to be 
unimportant. 

Typically, for any lattice computation one needs to 
evaluate the expectation value of an observable 0, 
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J _DC/cxp(-SG)e(m„) Dct M{m^) 
J DU cxp{~Sg) Dct M(m^) 



(1) 



where M is the Dirac matrix in x, colour, spin, flavour 
space for sea quarks of mass nis, Sg is the gluonic ac- 
tion, and the observable O may contain fermion prop- 
agators of mass m„. Sq ~ (^^trUpiaq/ g^, with go the 
bare coupling and Upiaq the product links along a pla- 
quette as shown in Figure [1] Amongst the many meth- 
ods of evaluation of eq.([T|), numerical simulations stand 
out due to the ability to achieve the goal of removing 
the lattice scaffolding, i.e., taking the continuum limit 
a ^ 0. Using the two- loop /3- function, it is easy to 
show that 



M-a=—{glho) 



■'5(l + 0(5o') 



(2) 



defines the way a mass scales Af -a on the lattice changes 
as the bare coupling 50 (= 1/6//?) is changed. Here 5o 
and bi are the universal coefficients of the /3-function. 
Typically, one needs larger and larger lattice sizes as 
a ^ in order to keep physical volume fixed. 

Numerically, the (O) is computed by averaging over 
a set of configurations {U^{xy\ which occur with proba- 
bility (X exp(— S'G)-Det M . Thus the main problem is to 
generate the ensembles of such configurations with the 
desired probability distribution. Complexity of evalu- 
ation of Det M has lead to various levels of approxi- 
mations in the process of generation of configurations: 
the quenched approximation consists of sea quark mass, 
nis = 00 limit whereas the full theory should have low 
sea quark masses: m„ — with a moderately heavy 
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strange quark. The computer time required to obtain 
results at the same precision increases as the sea quark 
mass is lowered. 

2.2. Some Results from Lattice QCD 



1.5 

1.0 

OeV 
0.5 

0.0 



Q5Q" 



o A 



O: Real world 

O: Lattice 

o; Fix parameter 



cf>33 

□nap 



OOK 



0>Tt 



□3 23 

□DlPI^ IP 



1~ Octet Dec. cc bb 



T r 



T 1 r 



■I . I 1 r 

il Average 

Hadronic Jets 



T r 



e'^e" rates 

1 Photo-production 

1 Fragmentation 
-o- 



1 Z widtli 



ep event sliape^ I 

rO-i- 



Polarized DI3 



Deep Inelastic Scattei;ing^(DIS) 

I ' T decays 
Spectroscopy (Lattice) 
Y decay 



J L 



J I li. 



J I L 



J L 



0.1 



0.12 



0.14 



Figure 3. Various determinations of a^. From [[5]. 



Figure 2. Comparison of experimental hadron spectra 
with lattice results [[4]. 

A variety of qualitative and quantitative results have 
been obtained using the lattice techniques. It will be 
both impractical and unnecessary to review all of them 
here. However, in order to appreciate the power of these 
techniques, we limit ourselves to providing a glimpse 
of them for the staggered fermions; similar, sometimes 
better in quality/precision, results have been obtained 
with the Wilson fermions as well. Figure [2] shows [ 
S] the results of the MILC and HPQCD coUaborations 
for the light as well as heavy hadrons obtained with 
light sea quarks. Using the pion and kaon masses to 
fix the scales of the corresponding quark masses, most 
other particle masses are found to be in good agree- 
ment with the experiment. Furthermore, the sponta- 
neous breaking of the chiral symmetry has been demon- 
strated by many groups since the early days of the lat- 
tice QCD, showing a non-vanishing chiral condensate, 
{■{pip) 0. Moreover, the goldstone nature of the pion 
has also been verified by checking that cx niu- Fig- 
ure [3] displays a comparison [ "5] of the lattice deter- 
mination of the strong coupling, as{Mz), with other 
perturbative determinations from experimental data. 



While these results verify that QCD is indeed the cor- 
rect theory of the strong interactions, and the lattice 
technique is the most reliable quantitative tool to ex- 
tract its non-perturbative properties, making new pre- 
dictions for the experiments is where the real challenges 
and excitement lies. It is very heartening to note that 
the decay constants of pseudo-scalar mesons containing 
a heavy quark were first obtained using lattice tech- 
niques: fD+ = 201±3±17MeVand /d^ = 249±3±16 
MeV [ [B]. These have since been measured experi- 
mentally to be fD+ = 223 ± 16 ± 7 MeV [ ^ and 
fo, = 283 ± 17 ± 14 MeV [[8], in excellent agreement 
with the lattice QCD predictions. 

2.3. Lattice QCD at Nonzero Temperature and 
Density 

Investigations of QCD under extreme conditions, 
such as high temperatures and/or densities, provide a 
solid platform for its most spectacular non-pcrturbative 
tests. Since the results from hadron spectroscopy fix 
the quark masses as well as the scale Aqcd , these tests 
are even completely free of any arbitrary parameters. 
Based on simple models, which build in the crucial 
properties of confinement or chiral symmetry break- 
ing and allow asymptotically for the free quark gluon 
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gas, one expects phase transitions to new phases such 
the Quark-Gluon Plasma or the colour superconduc- 
tors. As we shall see in the next section, the exper- 
imental possibilities of creating the required tempera- 
ture, and thus the new QGP phase, exist in the heavy 
ion collisions at high energies in BNL, New York and 
CERN, Geneva. Considering the scale of the entire ex- 
perimental enterprise, both in man-years invested and 
money spent, it seems absolutely necessary to have a 
better theoretical foundation for these results compared 
to merely relying on simple models. Fortunately, one 
can use the canonical Euclidean field theory formalism 
for equilibrium thermodynamics to look for the new 
phases, and the phase transitions in ab initio calcu- 
lations from the underlying field theory, i.e., QCD. In- 
deed, properties of the QGP phase can be predicted the- 
oretically using the lattice QCD approach, and tested 
in the experiments at BNL and CERN. As a first princi- 
ples based and parameter-free approach. Lattice QCD 
is an ideal reliable tool to establish the QCD phase dia- 
gram and the properties of its many phases. While most 
other basic features of the lattice formalism required for 
such an exercise remain the same as in section 2.1, a 
key difference for simulations at finite temperature is 
the need of an x Nt lattice with the spatial lat- 
tice size, Ng ^ Nt, the temporal lattice size for the 
thermodynamic limit V = N^a^ — > oo. The temper- 
ature T — 1/ {Nt ■ a) provides the scale to define the 
continuum limit : Fixing the transition temperature in 
physical (MeV) units and using cq. ([2]), the continuum 
limit is obtained by sending Nt oo. 

The lattice QCD approach has provided information 
on the transition temperature, the order of the phase 
transition, and the equation of state of QCD matter. 
One exploits the symmetries of the theory to construct 
order parameters which are then studied as a function 
of temperature to look for phase transitions, if any. 
QCD has two different symmetries in opposite limits of 
the quark mass ruq. For Nf flavours of massless quarks, 
QCD has SU{Nf) x SU{Nf) chiral symmetry while for 
niq — » oo, it has a global Z{3) symmetry. Such symme- 
tries usually imply zero expectation values for observ- 
ables which transform nontrivially under it unless the 
symmetry is broken spontaneously due to dynamical 
reasons and the vacuum transforms nontrivially under 
it. Lattice techniques enabled us to establish that the 
chiral symmetry is broken spontaneously at low tem- 
peratures, as indicated by its non-vanishing order pa- 
rameter, the chiral condensate {ipi') ^ 0. Its abrupt 
restoration to zero at high temperature will be a signal 
of a chiral symmetry restoring phase transition. Since 
the chiral condensate can be regarded as an effective 



mass of a quark, arising due to QCD interactions, the 
chiral transition can be interpreted as thermal effects 
'melting' this mass. Similarly, the global Z(3) symme- 
try breaking can be shown to be equivalent to a single 
quark having a finite free energy, i.e., the existence of 
a free quark. A nonzero expectation value for its or- 
der parameter, the Polyakov loop (i), is the a signal 
for deconfinement. Of course, in our world with two 
light and one moderately heavy flavours, neither sym- 
metry is exact but these order parameters may still act 
as beacons for transitions, depending on how mildly or 
strongly broken they are. 

2.4. Results from Lattice QCD at T ^ 0. 

The transition temperature can be determined by 
locating the point of discontinuity or sudden change in 
the order parameter as a function of the temperature 
(or other external parameter such as density). Since 
numerical results are necessarily obtained on finite lat- 
tices, there is an inevitable rounding which makes the 
determination of a little tricky. A lot of work has 
been done on this question in the statistical mechanics 
area and standard finite size scaling techniques exist to 
pin down Tc as well as the order of the transition. Since 
the early days, numerical simulations of lattice QCD 
have progressively tried to approach the real world of 
light quarks with vanishing effects from the lattice cut- 
off. The efforts began from the quenched approxima- 
tion, i.e., QCD without dynamical quarks, where the 
deconfinement order parameter (L) on small iVj-lattices 
was used to establish a first order deconfinement phase 
transition. Later QCD with three or more light dynam- 
ical quarks was also shown to have a first order chiral 
transition. Recent work on simulations for QCD with 
a realistic quark spectrum seems [ "9^ to rule out a first 
order chiral transition or a second order transition with 
the expected 0(4)-exponents, but suggests a rapid cross 
over. Determination of Tc, now the point of sharpest 
change, is even more tricky as a result. The current 
range for it can be summarized to be 170-190 MeV. A 
value on the lower end of the range was obtained [llOj by 
using larger A^j-lattices while a value at the upper end 
was obtained [ lllj using improved action but smaller 
Nt- There are other technical differences, such as the 
physical observable used to set the scale of lattice QCD, 
as well. Since the energy density is proportional to T'*, 
the current uncertainty in the value of translates to 
a ~ 60 % difference in the corresponding energy density 
estimates at Tc. In view of the tremendous impact it 
has on the requirements of heavy ion collision experi- 
ments, it is hoped that a narrowing of the range takes 
place as a result of future lattice QCD work. 
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Figure 4. Energy density and Pressure from lattice 
QCD. Taken from [|12j. 

Quantities of thermodynamic interest such as the en- 
ergy density, or the pressure or various quark number 
susceptibiHties can be obtained by using the canonical 
relations from statistical mechanics. Thus, 



dlnZ 



OT 



" ^^^^^ 



, etc. (3) 



Early results in the quenched QCD showed the exis- 
tence of a QGP phase which has energy density of about 
85 % of the corresponding ideal gas. The progress since 
then has been in employing large Nt and inclusion of 
light quark loops. Figure [4] displays recent results from 
such efforts. Obtained on two different lattice sizes, 
Nt — 4: and 6 with nearly realistic u, d and s masses, 
these results also exhibit similar kind of, ~ 15%, devi- 
ations from the ideal gas and do seem to hint towards 
the lattice cut-off effects to be small . The spatial vol- 
umes are perhaps not large enough to ensure that the 
thermodynamic limit is reached. However, this ques- 
tion is likely be addressed in near future soon. The 
results also suggest at most a continuous transition or 
even a rapid cross over; a strong first order phase tran- 
sition assumed / constructed in many phenomenological 
models seems clearly ruled out. This has implications 
for the hydrodynamical models used to analyse the ex- 
perimental data: possible mixed state of quark-gluon 
plasma and hadronic gas must be short lived, if at all 
it exists. 

From a theoretical perspective investigation of equa- 
tion of state offers hints of developing analytic or semi- 
analytic approaches. Thus conformal invariant theories 
are known to yield a variety of predictions for the ther- 



Figure 5. Entropy density s (in units of ideal gas en- 
tropy So) as a function of 't Hooft coupling. From [ 

[13]. 



modynamic quantities using the famous AdS-CFT cor- 
respondence. FigureOshows an attempt to confront the 
entropy density [ [13 for the quenched QCD in terms 
of the entropy of the ideal gas with the prediction of 

= 4 SYM [mj. The agreement is impressive, consid- 
ering the differences of the underlying theories. On the 
other hand, it is really in the stronger coupling region 
that it is not as good. Moreover, resummed weak cou- 
pling perturbation theory approaches seem to perform 
equally well at the lower couplings. Figure [5] shows the 
results [ T3 for the equation of state to highlight how 
conformal QCD really is. The ellipses denote 66% er- 
ror bounds on the measured EOS. The wedges piercing 
the ellipses have average slope c^, the speed of sound 
and the opening half-angle of these wedges indicate the 
error in . Conformal invariance is indeed violated sig- 
nificantly in the region close to the transition, with least 
violation at the same temperatures where in AdS-CFT 
prediction does well in Figure (5] 

Viscosities of the quark-gluon plasma, both the shear 
(77) and bulk (C), can also be determined using the 
lattice approach although unlike the equation of state 
these determinations need extra ansatze some of which 
are not universally accepted. Kubo's linear response 
theory lays down the framework to obtain such trans- 
port coefficients from certain equilibrium correlation 
functions. In particular, one obtains correlation func- 
tions of energy-momentum tensor using the lattice ap- 
proach above. These are, of course, defined at discrete 
Matsubara frequencies. Recall that the simulations at 
T ^ need lattices with i) periodic boundary condi- 
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Figure 6. Equation of State for (quenched) lattice 
QCD. Taken from [[13]. 



tions and ii) small Nt compared to Ng . The correlation 
function is thus defined at few discrete points only. One 
then continues it analytically to get the so-called re- 
tarded propagators in real time from which the the i] 
and ^ are obtained in the zero frequency limit. Fig- 
ure [7] shows the results [ |I5j in the quenched approxi- 
mation. Close to Tc, rather small values are obtained 
for the ratio of 77 to the entropy density s. These are 
seen to be consistent with the famous bound [T6^ from 
AdS-CFT. As shown in the Figure, perturbation the- 
ory suggests rather large values for this ratio. These 
results have since been refined [ I17j and made more 
precise but the general picture remains the same, as 
do the various theoretical uncertainties which plague 
these determinations. Larger lattices and inclusion of 
dynamical quarks will surely reduce some of these in 
near future. What is needed though for a more con- 
vincing demonstration of the fact the shear viscosity is 
indeed as small as hinted by the experimental data (see 
the next section) is a better control over the systematic 
errors in the analytic continuation. 

Analogous to the baryon number susceptibility, de- 
fined in eq. ([3]), various quark number susceptibilities 
can be defined by taking derivatives with the appro- 
priate chemical potential. These determine the fluc- 
tuations in the given conserved quantum number, say, 
strangeness. It has been argued [[T^ that under certain 
assumptions, testable experimentally, the strange sus- 
ceptibility can be related to the Wroblewski parameter 
As extracted from the data of heavy ion collisions. In- 
terestingly, lattice QCD computation in both quenched 
approximation and full QCD yield a As(Tc) ~ 0.4 — 0.5, 



Figure 7. Ratio of shear viscosity to entropy in 
(quenched) QCD vs. temperature. Taken from [il5j. 

whereas various experimental results [ [19] lead to a 
value 0.47 ± 0.04. Taking derivatives with two different 
chemical potentials in eq. ([3]), one obtains off-diagonal 
susceptibilities. These have the information on flavour 
correlations. Such a baryon-strangeness [[^D] or electric 
charge-strangeness [ [18] correlation has been proposed 
as a signature for identifying the nature of the high 
temperature phase as that of the quark-gluon phase. 
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Figure 8. Baryon-Strangeness and Electric charge- 
Strangeness correlation vs. temperature [I18j. 

Figure [5| shows the lattice results for QCD with 2 light 
dynamical quarks for both these correlations. They 
have been so normalized that a value unity, as seen in 
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most of the high temperature phase in Figure [51 char- 
acterises the existence of quark degrees of freedom with 
the appropriately fractional baryon number or charge. 
It has been shown that the correlation in the low tem- 
perature phase are consistent with the hadronic degrees 
of freedom. Indeed, any lack of the expected transition 
should lead to much milder temperature dependence 
as well as a value different from unity for these cor- 
relation functions. Being ratios of the quark number 
susceptibilities, these correlations are robust, both the- 
oretically and experimentally. Systematic errors due to 
lattice cut-off or dynamical quark masses are therefore 
very small as are the systematic errors from experimen- 
tal sources. 
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Figure 9. Debye radii for charmonia vs. temperature [ 

ED- 

Debye screening of coloured heavy quarks in the de- 
confined phase had long been recognised [123] as a pos- 
sible signal of formation of quark-gluon plasma, de- 
tectable in the suppression of heavy quarkonia in the 
heavy ion collisions. In view of the impressive data 
from CERN at lower SPS energies, and the expec- 
tations from the upcoming LHC experiments, a crit- 
ical assessment of the original theoretical argument 
seems prudent. Lattice QCD has contributed hand- 
somely in finite temperature investigations of both the 
heavy quark-antiquark potential, which can be used 
in the Schrodinger equation to look for the melting of 
heavy quarkonia, and directly in the spectral function 
at finite temperature. Figure [9] displays the results [ 
I21j for the screening radii estimated from the inverse 
non-perturbative Debye mass mu in quenched (open 
squares) and full (filled squares) QCD. For r < r„ied, 



Figure 10. Spectral function of rjc and J/ip- From [ 

Eg. 

the medium effects are suppressed, leading to the same 
heavy quark potential as at T = 0. The horizontal 
lines correspond to the mean squared charge radii of 
J/ '4' 7 Xc and ipf charmonia, and are thus the averaged 
separations r entering the effective potential in poten- 
tial model calculations. Figure [9] therefore suggests 
that the Xc and ■(/;/ states would melt just above the 
transition while J/^ may need higher temperatures to 
be so affected. Direct spectral function calculations 
[ I22j provide a strong support for such a qualitative 
picture. Such computations have been made feasible 
by the recognition of the maximum entropy method 
(MEM) technique as a tool to extract spectral func- 
tions from the temporal correlators computed on the 
Euclidean lattice. However, as in the case of shear vis- 
cosity above, the data for such temporal correlators are 
sparse, making the extraction more of an art. Neverthe- 
less, large lattices, 48'^ x 12 to 64^^ x 24 have been used in 
this case to avoid such criticisms. Figure [TOl shows typ- 
ical results for the J/^p and rjc mesons in the quenched 
approximation. The vertical error bars denote the pos- 
sible uncertainties on the area under the peak as defined 
by the horizontal error bar. The peaks in both spec- 
tral functions appear to persist up to 2.25 Tc, i.e., have 
nonzero area within the computed error-band, and are 
gone by STc unlike the Xc which has no peak already by 
1.1 Tc- Further technical improvements, such as the in- 
clusion of light dynamical quarks, are clearly desirable. 
Another important issue is that of the huge widths of 
the peak compared to their known zero temperature 
values. If real, they could hint at rather loosely bound 
states which could be dissociated by thermal scatter- 
ings. 
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2.5. QCD Phase Diagram 

The quark-gluon plasma phase and the correspond- 
ing quark-hadron transition which we discussed so far 
is a special case of the conditions that could be created 
in the heavy ion collisions. Indeed, the lattice QCD 
thermodynamics that we considered was for the case 
of zero net baryon density and an almost baryon-free 
region can be produced in the heavy ion collisions in 
the so-called central rapidity region, as wc explain in 
the next section. It also pervaded our Universe a few 
microseconds after the Big Bang. In general, of course, 
one should expect hot regions with some baryon num- 
ber since the colliding nuclei themselves carry substan- 
tial baryon number. Massive stars could also have re- 
gions of huge baryon densities in the core which could 
even be at rather low temperatures. It is natural to 
ask what these generalized extreme conditions lead us 
to. One could have new phases, and different natures 
of phase transitions which may even have astrophysical 
consequences. The vast research area of QCD phase 
diagram in the plane of temperature T and the bary- 
onic chemical potential fj,B deals with these and several 
other interesting issues. While the current theoretical 
expectations suggest such physics at nontrivial baryon 
densities to be better accessible to the colliders at lower 
energies, such at the RHIC in New York or the forth- 
coming FAIR facility at GSI, Darmstadt, we feel that 
the physics may be interesting in its own right to be in- 
cluded in this article dedicated to LHC; with some luck 
LHC experiments may have important contributions to 
this area as well. 

Using simple effective QCD models, such as the 
Nambu-Jana Lasinio model at finite temperature and 
densities [ [24], several speculations have been made 
about how the QCD phase diagram in the T-fiB plane 
should be. At asymptotically high densities, one ex- 
pects quarks to be effectively free, and therefore to ex- 
hibit various colour superconducting phases [ [5S] . In 
the limit of large number of colours Nc for quarks, it has 
also been argued that a "quarkyonic" phase may exist 
[ [26] at low enough temperatures. A crucial question, 
especially in the context of either the massive stars, or 
heavy ion collisions, is the quantitative reliability of the 
predicted regions in the T-fiB space. Alternatively, it 
is unclear how low can the asymptotic predictions be 
trusted. Nevertheless, most model considerations seem 
to converge [[25] on the idea of the existence of a critical 
point in the T-^b plane for the realistic case of 2 light 
flavours (m„ — rud) of dynamical quarks with a moder- 
ately heavy strange quark. Establishing it theoretically 
and/or experimentally would have huge profound con- 
sequences in our (non-perturbative) understanding of 



QCD. 

Extending the lattice approach to the case of QCD at 
finite density has turned out to be a challenging task at 
both conceptual and computational level. In principle, 
it really is straightforward. One just has to add a term 
(J-bNb = ^B'4'70^ term to the fermionic part of the 
action, hence the Dirac matrix M, in eq.([T]). In order 
to eliminate certain spurious divergences, even in the 
free case, some care is needed [ 27 and the naive form 
above has to be modified. A big conceptual block has, 
however, turned up in form of our inability to define 
exact chiral invariance in the presence of the chemi- 
cal potential [ 28J : both the Overlap and the Domain 
Wall fermions lose their exact chiral invariance for any 
nonzero fi. The staggered fermions do preserve the chi- 
ral invariance for nonzero /i. Furthermore, they are sim- 
pler to handle numerically. Again most of the numerical 
work has therefore employed the staggered fermions, al- 
though they are plagued with the difficulties of precise 
definition of flavour and spin as mentioned earlier. In- 
deed, the existence of the critical point depends [ [^ 
crucially on how many flavours of light quarks the the- 
ory has. Proceeding none the less with the staggered 
quarks, another tough problem arises in form of the fact 
that the Det M(/i ^ 0) in eq. H]) is complex whereas 
the numerical methods of evaluation, employed to ob- 
tain the results in the sections above, work only if the 
determinant is positive definite. This is akin to the sign 
problem well known to the statistical physicists and is 
largely unsolved in its full generality. 

1.1 , , , , , , 



1 




Figure 11. QCD Phase diagram for 2 light fiavours of 
quarks. The circles [[29l[3T] and the square [[32] denote 
the location of the critical point on lattices with 1 /4T 
and 1/6T cut-offs respectively. Taken from ["iTr, where 
more details can be found. 
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A bold breakthrough was achieved [ [53] by applying 
the method of re-weighting in the vicinity of the finite 
temperature transition at /i = 0. A flurry of activ- 
ity saw many new methods emerge [ 30^, such as ana- 
lytic continuation of computations at imaginary chem- 
ical potential and Taylor series expansions of the free 
energy. These have been employed to get a glimpse 
of whether a critical point does exist, and if yes, what 
its location may be. The field is really in its infancy 
and unfortunately at present no consensus amongst the 
results obtained so far has emerged. Figure [11] ex- 
hibits the results obtained for the critical point for the 
case of two flavours of light quarks with a pion mass 
^tt/itT'p — 0-31 ± 0.01, compared to 0.18 in the real 
world. The results [ [221 [31] denoted by circles in the 
Figure [TT] are for a lattice cut-off a = 1 /4T whereas 
the square [ [32] denotes the first attempt towards the 
continuum limit by lowering a to 1/6T. Large finite 
volumes have been observed. The shift in the location 
of the open circle in the Figure [TT] was shown ['3T] to 
be due to the use of a 10 times larger volume than the 
open circle In order to be brief, we prefer to close 

this section by noting that different results have been 
claimed in the literature for larger pion masses and for 
a different number of flavours. It is hoped that a clear 
and solid picture will emerge in the near future. 

3. Relativistic Heavy-Ion Collisions 

At energies of a few GeV/N to a few lO's of GeV/N, 
colliding nuclei tend to stop each other thereby form- 
ing a dense, baryon-rich matter. At higher ener- 
gies, they nearly pass through each other forming a 
dense, nearly baryon-number-free matter in the mid- 
rapidity region. This is evident in the shapes of ra- 
pidity distributions {dN/dy vs y) of the net proton 
(i.e., proton— antiproton) production observed at var- 
ious beam energies. This apparent transparency of nu- 
clear matter at ultra-relativistic energies can be under- 
stood in the space-time picture of the collision, pro- 
posed by Bjorken [[55] [M]. 

3.1. Bjorken Picture 

Consider, for simplicity, a central (i.e., head-on or 
zero impact parameter) collision of two identical spher- 
ical nuclei in their CM frame. Coordinate axes are cho- 
sen such that the two nuclei approach each other along 
the 2;-axis and collide at the origin at time t = Q. Deep 
inelastic scattering experiments have revealed the par- 
ton structure of hadrons: In the proton, e.g., the va- 
lence quark distributions a:u„(a;), xdv{x) peak around 
a; ~ 0.2 and vanish as x ^ 0/1. (x is the Bjorken scal- 



ing variable.) The gluon and sea quark distributions, 
xg{x), xUs{x), xds(x), on the other hand, shoot up as 
X 0. These numerous low-momentum partons are 
called wee partons. As a result of the Lorentz contrac- 
tion, the longitudinal (i.e., parallel to the beam axis) 
spread of the valence quark wave function is reduced to 
^ 2i?/7 where R is the nuclear radius and 7 its Lorentz 
factor. However, no matter how high the beam energy 
(or 7), the incoming nuclei always have in them wee 
partons with typical momenta p ~ ^qcd, and hence 
longitudinal spread ~ 1 fm [ '33^ . The wee partons pre- 
vent the nucleus from shrinking below ~ 1 fm in the 
z-direction. If 2R/"/ < 1 fm, they play an important 
role in the collision dynamics. 

As a result of the collision of two nuclei, or rather 
two clouds of wee partons, a highly excited matter 
with a large number of virtual quanta is created in 
the mid-rapidity region. (In the modern parlance one 
talks about coherent "glasma" formed by a collision 
of two sheets of "colour glass condensates (CGC)" [ 
I35j.) Hereinafter we discuss only the mid-rapidity re- 
gion. The virtual quanta need a finite time (t^cc) to 
decohere and turn into real quarks and gluons. Here 
Tdoc refers to the rest frame of an individual parton. In 
the overall CM frame, the relevant time is "/t^cc due to 
the time dilation, 7 being the Lorentz factor of the par- 
ton. It is now clear that "slow" partons decohere earlier 
and hence near the origin, than the "fast" ones which 
emerge later at points farther away from the origin. 
(This is known as the inside-outside cascade.) In other 
words, the large- a; part of each nuclear wave function 
continues to move along its light-cone trajectory leav- 
ing the small- a; part behind. Thus, in the limit of high 
beam energy, the time dilation effect causes the near 
transparency of nuclei, referred to earlier. 

Figure [T2I shows this schematically in 1-1-1 dimension 
for simplicity. The curves are hyperbolas of constant 
proper time r = y/W^^z^. All points on a given hyper- 
bola are at the same stage of evolution. In particular, 
let the hyperbola labelled '1' refer to r = Tdoc- Parton 
at z undergoes decoherence at time t = \/t^^^ + z^. 
The larger the z, the larger the time t and higher the 
parton velocity — z/t [ 34]. 

If the partons thus formed interact amongst them- 
selves a multiple number of times, the system ap- 
proaches local thermal equilibrium. Thermalization 
time Tth (> Tdcc) is estimated to be of the order of 
1 fm. 

Figure [T^ indicates a possible scenario. 1,...,5 are 
the hyperbolas with proper times ri, ...,T5. 
t — — z : the instant of collision 
< r < Ti : formation of quark-gluon matter 
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Ti < T < T2 : (local) equilibration of quark-gluon 
matter, i.e., formation of QGP 

T2 < r < T3 : hydrodynamic evolution of QGP 
(partonic EOS) 

T = T3 : hadronization 

T3 < T < T4 : hydrodynamic evolution (hadronic EOS) 
T4 < T < T5 : transport theoretic evolution of hadrons 
T — : freezeout 

T > T5 : free-streaming to detectors 
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Figure 12. Space-time picture of an ultra-relativistic 
nucleus-nucleus collision in 1 -|- 1 D for simplicity 



The above is a rather simple-minded picture: in real- 
ity, there are no such "water-tight compartments" . The 
framework of hydrodynamics is applicable, if at all, only 
when the system is at or near (local) thermal equilib- 
rium. If the matter formed in ultrarelativistic heavy- ion 
collisions is fully thermalized, one may use the frame- 
work of relativistic ideal fluid dynamics to study its 
evolution. If it is only partially thermalized, one could 
use relativistic dissipative fluid dynamics. In any case, 
the covariant transport theory provides a more general 
framework for this purpose. 

Bjorken [ [34j presented the following formula to es- 
timate the energy density attained in the mid-rapidity 
region: 

1 AEt 

^0 = ^2 1- ' (4) 

where R is the nuclear radius, tj ~ 1 fm/c is the for- 
mation time of QGP, and Et is the transverse energy. 



It is clear that even if QGP is formed, its lifetime will 
be of the order of a few fm/c or 0(10"^'^) seconds, and 
what experimentalists detect in their detectors are not 
quarks or gluons, but the standard hadrons, leptons, 
photons, etc. It is a highly nontrivial task to deduce the 
formation of QGP from the properties of the detected 
particles. This is analogous to the situation in cosmol- 
ogy where one tries to deduce the information on the 
early epochs after the Big Bang by studying the cosmic 
microwave background radiation and its anisotropy. 

Actually the analogy between the Big Bang and the 
"Little Bang" is quite striking. In both the cases the 
initial conditions are not accurately known, but there 
are plausible scenarios. In the former case, there is 
inflation occurring at ~ 10"'^^ sec, with the inflaton 
energy converting into matter and radiation, leading to 
a thermal era. In the latter case, one talks about a 
highly excited but coherent glasma converting, on the 
time scale of ^ 10^^"* sec, into quarks and gluons which 
may thermalize to form QGP. In both the cases the 
"fireball" expands, cools, and undergoes one or more 
(phase) transitions. Decoupling or freezeout follows — 
of photons in the former case and of hadrons in the lat- 
ter. The unknown initial conditions are parameterized 
and one tries to learn about them by working one's way 
backwards, starting with the detected particles. As we 
shall see shortly, the anisotropy of the detected parti- 
cles plays a crucial role in the diagnostics of the Little 
Bang too. 

Definition: The STAR collaboration at RHIC has 
defined the QGP as "a (locally) thermally equilibrated 
state of matter in which quarks and gluons are decon- 
fined from hadrons, so that colour degrees of freedom 
become manifest over nuclear, rather than merely nu- 
cleonic, volumes" [|36]. The two essential ingredients 
of this definition are (a) local equilibration of matter, 
and (b) deconfinement of colour over nuclear volumes. 
Recent claims of the discovery of QGP at RHIC [ [37] 
were based on two observations which, for the first time, 
provided a good evidence that each of these two require- 
ments has been fulfilled. We discuss them one by one 
in the next two subsections (3.2, 3.3). That will be 
followed by brief descriptions of a few other signals of 
QGP in subsections 3.4, 3.5. 

3.2. Anisotropic Flow 

Consider now a non-central (or non-zero impact pa- 
rameter) collision of two identical (spherical) nuclei 
travelling in opposite directions. Choose x, y axes as 
shown in Fig. 1131 The collision or beam axis is per- 
pendicular to the plane of the figure. Length of the 
line AB connecting the centres of the two nuclei is the 
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impact parameter b. Plane xy is the azimuthal or trans- 
verse plane. Plane xz is the reaction plane. It is de- 
termined by the impact parameter vector b and the 
collision axis. (Obviously the reaction plane cannot 
be defined for a central collision.) — tan~^ (py/p^) 
is the azimuthal angle of an outgoing particle. The 
almond-shaped shaded area is the overlap zone. In a 
real experiment, Fig. I14i the x, y axes need not coin- 
cide with the lab-fixed X, Y axes. Indeed the reaction 
plane subtends an arbitrary angle (f>R with the X axis. 
(f)R varies from event to event. It is a priori unknown 
and special experimental techniques are needed for its 
determination. 




Figure 13. Non-central collision 




^ X 



Figure 14. Non-central collision. XY are lab-fixed axes. 



The triple differential invariant distribution of parti- 
cles emitted in the final state of a heavy-ion collision 
is a periodic even function of (p, and can be Fourier 
decomposed as 



E- 



(Pp 



PTdpTdyd4i 
d^N 1 
pxdprdy 27r 



1 



2vn cos(n0) 



where y is the rapidity and (j) is measured with re- 
spect to the reaction plane. The leading term in the 
square brackets in the above expression represents the 
azimuthally symmetric radial flow, vi is called the di- 
rected flow and V2 the elliptic flow. w„ = (cos(n0)) is 
actually a function of pt and y. Here the average is 
taken with a weight equal to the triple differential dis- 
tribution of particles in the (pT, y) bin under consider- 
ation. V2 can also be written as {{p^ — pI)/{Px ~^Py))- 
For a central collision the distribution is azimuthally 
isotropic and hence w„ = for n — 1,2,.... In other 
words, only the radial flow survives. 

Measurement of the radial flow: Radial flow gives a 
radially outward kick to the emerging hadrons thereby 
depleting the low-px population and making their pT 
spectra flatter. The heavier the hadron, the stronger 
the momentum kick it receives. By measuring the 
slopes of the px spectra of various hadrons, the radial 
flow velocity can be extracted. At RHIC it turns out 
to be a sizeable fraction (~ 50%) of the speed of light. 
Thus the flow is compressible. 

Measurement of the anisotropic flow w„: There are 
several methods, (a) The most obvious one is based 
on the deflnition w„ = {cosn{(j)— (f>R)) where both (j) 
and are measured with respect to a lab-fixed frame 
of reference. This, however, requires the knowledge of 
which varies from event to event and is not easy to 
determine, (b) Two-particle correlation method: This 
gives v'^ = (cosn((/)i — 02)), where 0i and <p2 are az- 
imuthal angles of two outgoing particles. This method 
has an advantage that the reaction plane need not be 
known. However, w„ is determined only up to the sign. 
There are several other methods such as the cumulant 
method [ [38] , mixed-harmonic method [ [39] , Lee- Yang 
zeroes method [[40], etc. For a recent review, see [|4T]. 

Importance of the anisotropic flow v„: Consider a 
non-central collision, Fig. [131 Thus the initial state is 
characterized by a spatial anisotropy in the azimuthal 
plane. Consider particles in the almond-shaped overlap 
zone. Their initial momenta are predominantly longi- 
tudinal. Transverse momenta, if any, are distributed 
isotropically. Hence ?;„ (initial) = 0. Now if these 
particles do not interact with each other, the final 
(azimuthal) distribution too will be isotropic. Hence 
u„ (final) = 0. 

On the other hand, if these particles interact with 
each other a multiple number of times, then the (local) 
thermal equilibrium is likely to be reached. Once that 
happens, the system can be described in terms of ther- 
modynamic quantities such as temperature, pressure, 
etc. The spatial anisotropy of the almond-shaped over- 
lap zone ensures anisotropic pressure gradients in the 
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transverse plane. This leads to a final state character- 
ized by a momentum anisotropy in the PxPy plane or 
equivalentljQ to an anisotropic distribution of particles 
in the transverse {xy) plane, and hence a nonvanishing 
Vn- Thus Vn is a measure of the degree of thermaliza- 
tion of the matter produced in a noncentral heavy-ion 
collision. 

To sum up, if either of the two ingredients, namely 
initial spatial anisotropy and adequate rescatterings, is 
missing, there is no anisotropic flow (w„). 

Sensitivity of w„ to properties of matter at early times 
fni/c): We saw above that the spatial anisotropy of 
the initial state (together with multiple rescatterings) 
leads to more matter being transported in the direc- 
tions of the steepest pressure gradients, and thus to a 
non-zero Vn- That in turn results in the reduction in 
spatial anisotropy ("self-quenching"). In other words, 
expansion of the source gradually diminishes its spa- 
tial anisotropy. Thus w„ builds up early (i.e., when the 
spatial anisotropy is significant) and tends to saturate 
as the spatial anisotropy continues to decrease. (This 
is unlike the radial flow which continues to grow until 
freeze-out and is sensitive to early- as well as late-time 
history of the matter). Thus u„ is a measure of the 
degree of thermalization of the matter produced early 
in the collision. In other words, u„ is a signature of 
pressure at early times. 

Hydrodynamic calculations of v„ involve the equa- 
tion of state of QGP. Thus one hopes to learn about the 
material properties of the medium, such as the speed of 
sound, sheer and bulk viscosities, relaxation times, etc. 

Flow may also be affected by the dynamics of the 
hadronic phase. Study of the flow would provide con- 
straints on the properties of hadronic matter too. (It is 
expected that at LHC, the relative contribution of the 
QGP phase to Vn would be larger than that at SPS and 
RHIC. This would reduce the effect of the uncertainties 
in the hadronic phase). 

It should, however, be kept in mind that the ini- 
tial conditions for the hydrodynamic evolution are not 
known with certainty. Hence the task of unravelling the 
properties of medium is not as easy as it may appear. 

Figure [T5] shows the impressive agreement between 
RHIC data on V2{pt) and ideal hydro calculations for 
Pt up to ~ 1.5 GeV/c. In particular note the mass or- 
dering: the heavier the hadron, the smaller the V2{pt)- 
This can be understood heuristically as follows. 

Mass ordering of V2{pt)- Recall that the radial flow 
depletes the population of low-p^ hadrons (by shift- 
ing them to larger values of pt)- This effect is more 

-•-Since (f) = tan~^{py/px). 



pronounced for larger flow velocities and for heavier 
hadrons. Suppose V2 is positive as at RHIC, which 
means more hadrons emerge in-plane (x-direction) than 
out-of-plane (y-dircction). Now due to higher pres- 
sure gradients in the a;-direction, hadrons which emerge 
in-plane experience a larger flow velocity than those 
which emerge out-of-plane. So the depletion is greater 
for the hadrons emerging in-plane than out-of-plane. 
This tends to reduce the anisotropy and hence V2 
of all hadron species. For a heavier hadron species 
this reduction is more pronounced. The net result is 

liqht hadron/ \ ^ heavy hadron/ \ -\ j j • 

V2 [Pt) > V2 \Pt)- Mass-ordermg 

signifies a common radial velocity field. 

Hydrodynamic model calculations predicted mass or- 
dering of V2{pt)- The broad agreement between the 
RHIC data and the predictions of ideal hydro (Fig. [T5|) 
led to the claims of thermalization of matter and dis- 
covery of a perfect fiuid — more perfect than any seen 
before. 

In order to claim the discovery of a new state of mat- 
ter, namely quark-gluon plasma, one needs to demon- 
strate unambiguously that (local) equilibrium is at- 
tained. There are indications that the equilibrium at- 
tained at RHIC is incomplete [ i42j. 



Vs>,., = 200GeV ^^^Au + ^^^Au at RHIC 



_ STAR PHENIX 
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Transverse momentum p, (GeV/c) 



Figure 15. Minimum-bias data. Curves represent 
ideal hydro results with a first-order QGP-hadron phase 
transition. Figure taken from [ 143] . 
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3.2.1. Constituent Quark Scaling 

For pt ^ 2 GeV/c, ideal hydro results are in gross 
disagreement with the V2{pt) data: calculated V2{pt) 
continues to rise with px^ while the data tend to satu- 
rate and the mass ordering is reversed. In the interme- 
diate momentum range (2 GeV/c ^ pt ^ 5 GeV/c), 
it is observed that the V2/nq vs pr/nq (or KET/nq) 
data fall on a nearly universal curve; see Fig. [161 Here 
Uq is the number of constituent quarks and KEt is the 
transverse kinetic energy. This is called the constituent 
quark scaling. It shows that the flow is developed at 
the quark level, and that the hadronization occurs by 
quark recombination. 
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Figure 16. Left: Note the two distinct branches. Right: 
Universal curve. Figure taken from [I44j. 



3.3. Jet Quenching 

A variety of signatures of quark-gluon plasma have 
been proposed. Some of the more popular ones are 
excess strangeness production, thermal dileptons and 
photons, jet quenching, J/i/i-suppression and event-by- 
event fluctuations. A common theme underlying all of 
these is the idea of exploiting the consequences of those 
properties of QGP which distinguish it from alterna- 
tives like a hot hadron gas. Since QGP is expected to 
form and exist predominantly in the early phase of the 
collision, the so-called hard probes are potentially the 
cleaner direct probes of this early phase. It is experi- 
mentally known that rare but highly energetic scatter- 
ings produce jets of particles : g + g g + g, where en- 
ergetic gluons from the colliding hadrons produce two 
gluons at large transverse momenta, which fragment 
and emerge as jets of showering particles. Their typical 
production time scale is t ^ where Q = pt, the 

transverse momentum of the jet, is the hard scale of 



production. Thus jets at large transverse momenta are 
produced very early and by traversing through the pro- 
duced medium carry its memory while emerging out. 
Quark-Gluon Plasma, or any medium in general, in- 
teracts with the jet, causing it to lose energy. This 
phenomenon goes by the name of jet quenching. 

Using the well-known factorization property of per- 
turbative QCD ['45 , which allows a separation between 
the hard and soft scales, a typical cross section at hard 
scale, say that of hadron h at large transverse momenta 
in the process A+B h+X, can be symbolically writ- 
ten as 



fAixi,Q^)(g 



(g)D,^h{z,Q^) . (5) 

Here /a, /b are parton distribution hmctions of the col- 
liding hadrons A and B at scale Q^, a{xi,X2, Q^) is the 
elementary pQCD cross section for partons of momen- 
tum fractions xi and X2 to produce a parton i with the 
hard scale Q = px for jet production, and Di^hiz, Q^) 
is its fragmentation function to hadron h with momen- 
tum fraction z. Various convolution integrations are 
denoted symbolically by (8". Clearly, there are many 
more details which are not spelt out here for brevity, 
such as the kinematic integration region or the summa- 
tion over all allowed many parton level processes, such 
quark-quark or gluon-quark etc. These can be found in 
textbooks [45!. 
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Figure 17. Comparison of the various dihadron angular 
correlations. Taken from [147). 



In presence of a medium, of hot hadron gas or quark- 
gluon plasma, the function D above will get modified 
by the interactions with medium. The medium pro- 
vides scattering centers for the fast moving seed particle 
of the jet which typically impart a transverse momen- 
tum kick to it. The medium induced transverse mo- 
mentum squared per unit path length, q, characterizes 
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Figure 18. Comparison of the dihadron azimuthal 
correlations as a function of the associated pr for 8 
< pp^^ < 15 GeV. Taken from [|48]. 



the quenching weight function P{AE) [ [46| which is 
the probabihty that a hard parton loses an additional 
energy AE due to its interactions with the medium. 
In hot matter with a temperature of about T = 250 
MeV, a perturbative estimate [ [49] for q is about 0.5 
GeV^/fm. It is typically a lot smaller in the cold nu- 
clear matter. In terms of the quenching weight, one 
can write down [[IB] a medium modified fragmentation 
function for a jet passing through a medium as 



l-e 



l-e' 



(6) 



For a heavy quarkoniuni like J /tp, the analogue of Z?, 
is the wave function of a heavy quark-antiquark pair 
(cc), and it will be presumably flatter in a hot medium, 
corresponding to "its melting" . 

RHIC experiments have cleverly exploited their capa- 
bilities to perform tests which have an on-off nature and 
are therefore rather convincing about the qualitative 
existence of the jet quenching phenomenon in heavy ion 
collisions. In the case of the elementary g + g ~^ g + g 
hard process, one expects back-to-back jets, i.e, a well- 
determined azimuthal correlation between the fast par- 
ticles. As jets are hard to identify in the complex multi- 
particle environment at RHIC, the STAR collaboration 
constructed the angular correlation of hadrons, using a 
high transverse momentum p^*^^ particle as the trigger, 



and studying the azimuthal distributions of the associ- 
ated particles (p!^"^"'^ < p^*^^). Figure flTl compares 
the results for gold-gold central collisions, where one 
expects formation of a hot medium, with the proton- 
proton or deuterium-gold collisions, where one expects 
to have turned off the medium effects. The expected 
correlation, signalling a lack of any quenching/medium, 
is clearly visible in the two peaks separated by 180° for 
the d-Au and pp collisions. Remarkably the gold-gold 
central collision data show only the peak at zero degree 
or the near-side. A hint of the creation of some medium 
is given by the vanishing of the away-side jet, at 180° 
degrees, which appears to have been fully quenched by 
the medium. For high enough trigger px, one can do 
the same comparison as a function of range of the as- 
sociated pt- Clearly, as the increases, one ought 
to see the away-side re-emerge. This is beautifully seen 
in the Figure [181 R shows the azimuthal correlations 
for 8 < p^*^^ < 15 GeV for d-Au, and Au-Au collisions 
in two centrality bins, with the data for most central 
collisions displayed in the last column. The pT of the 
associated particle is restricted to ranges marked on 
the right side, and increases as one goes from top to 
the bottom. All panels show comparable strengths for 
the near-side peak. As the p^^"^"'^ grows above 6 GeV, 
the away-side peaks in all the three systems also show 
comparable strengths whereas for lower p^'*''°'^ ranges 
one has diminishing away-side peaks, characteristic of 
jet-quenching. The same phenomena can also be stud- 
ied by varying the p^'^^ and the away-side peak is seen 
clearly to emerge as p^'^^ increases. 

A more quantitative investigation of the jet quench- 
ing phenomena needs to extract the transport coefh- 
cient q, and establish the presence of the hot matter 
by comparing it with the corresponding theoretical es- 
timates, directly from QCD. Many such attempts have 
been made. Recently, the PHENIX experiment [ [5T] 
reported their measurement of neutral pion production 
in Au-Au collisions at 200 GeV at the RHIC collider in 
BNL. They define the now-famous nuclear suppression 
factor Raa as the weighted ratio of the nuclear dif- 
ferential distribution in rapidity y and transverse mo- 
mentum Pt and their own earlier measurements for the 
same quantity in proton-proton. 



1/N^^tdN/dydpT 
{TAB)d(Tpp/dydpT 



(7) 



where further details of determinations of various fac- 
tors above are given in [ [^ . Their results for Raa are 
displayed in Figure [THl While the first panel shows the 
results for their entire data set, the other panels exhibit 
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Figure 19. Nuclear modification factor, Raa, for neu- 
tral pions as a function of transverse momentum for 
different centralities. Taken from [I51|. 



data for increasing peripherality of the collisions (indi- 
cated by the increasing range of the percentage label 
of each panel) , or decreasing centrality. The error bars 
indicate the statistical errors, whereas various system- 
atic errors are shown by the boxes. Note that if the 
nucleus-nucleus collisions were merely scaled proton- 
proton ones, one expects Raa = 1- What the data 
in Figure [19] indicate, however, is a five-fold suppres- 
sion that is essentially constant for 5< pr <20 GeV for 
the most central bin of 0-10 %. The qualitative pattern 
is the same in all centralities, although the magnitude 
of suppression comes down. The highest centrality bin 
was used to determine the transport coefficient in the 
the parton quenching model [^SOj to obtain q = 13.2^^ \ 
GeV^/fm. Typically, fits with varying model assump- 
tions do tend to yield a g of 5-15 GeV^/fm. This order 
of magnitude or so higher value of the transport coef- 
ficient compared to the expectations from perturbative 
QCD, ~ 0.5, as mentioned above is an unresolved puz- 
zle. Nevertheless, the value hints at a hot medium, 
presumably even stronger interacting than the pQCD 
picture, as the cold matter expectations for q are even 
more in disagreement with the experimental determina- 
tion. Clearly a lot more needs to be understood from 
the data by further delving into the detail predictions 
of the models and confronting them with data, as [ I51| 
attempts to do, in order to establish the nature of hot 



medium produced as that of quark-gluon plasma. 

Having discussed the two main observations, 
anisotropic flow and jet quenching, which lend support 
to the claims of discovery of QGP at RHIC, we now 
discuss some corroborative evidences which strengthen 
these claims. There are also surprises in the RHIC data 
when compared with the expectations from the earlier 
lower energy heavy ion collisions at SPS in CERN. We 
discuss some with the aim to prepare ourselves for the 
expectation at yet higher energy in LHC. 

3.4. Anomalous J/ip Suppression 

Amongst the many signatures proposed to look for 
QGP experimentally, the idea of JZ-^-suppression has 
attracted the most attention as the likely "gold-plated" 
signal. Soon after the pioneering work of Matsui and 
Satz [|23', arguing that i) as a hard QCD process, the 
heavy charm pair production takes place very early, ii) 
the Debye screening of the QGP prevents formation of 
a J/ip state in heavy ion collisions, and iii) the low 
temperatures at the hadronization do not permit pro- 
duction of charm- anticharm pair kinematically, it was 
further proposed that the suppression pattern ought to 
have a characteristic [ [52] transverse momentum de- 
pendence. Recognising that the gluon and quark dis- 
tribution functions depend on the atomic number A, 
known by the famous EMC-effect, it was shown in a 
perturbative QCD calculation that the suppression sig- 
nal [[53] itself as well as its pr-dependence [ 54 can be 
mimicked by the mundane nuclear shadowing. Thus it 
became clear since the early days that a detailed quan- 
titative analysis is necessary to disentangle the effects of 
the Debye screening in QGP. It has since been recog- 
nised that other effects, notably the absorption [ [55] 
of the produced J/ip in the nucleus, causes suppres- 
sion of J/tp in all pA and A_B-collisions. Thus one has 
to first account for this expected or normal suppres- 
sion and then look for additional or anomalous J/ijj- 
suppression as the possible signal of QGP. Considering 
the general wisdom that J/V>production can be com- 
puted in pQCD, it ought to be a straightforward task to 
compute this normal suppression. Unfortunately, it is 
not so. One reason is that the gluon distribution func- 
tion, and the nuclear shadowing effects, are not well 
known. Another, perhaps much more important rea- 
son, is that the hadroproduction of J/tp needs to tackle 
the vexing issue of its formation from the perturbatively 
produced charm- anticharm pair. One usually depends 
[ 156] on models, such as the colour evaporation or the 
color octet model, hoping that the effective theory de- 
scriptions are valid. It turns out to be true for large 
Pt charmonium production but not for the total cross 
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Figure 20. J/'0-suppression in Pb-Pb collisions at SPS 
as a function of transverse energy Et- Figure taken 
from [[57]. 



The preferred phenomenological method [ I55j has 
been to parametrise the ratio of J/ip-cioss sections, to- 
tal or appropriate differential cross sections in its trans- 
verse momentum px , or forward momentum fraction x p 
etc., in pA and pp collisions at the same colliding en- 
ergy, y^, as exp(— crabs(^/V')poi), where L is the mean 
length of the trajectory of the produced cc pair in nu- 
clear matter and po is the nuclear density. The parame- 
ter, a^hsiJ/'ip), is obtained by fitting the data. Defining 
a mean free path A — l/usx^siJ/''P)Poi one then extends 
this idea to the heavy-ion collisions to define the normal 
or expected J/ip suppression due to the traversing of the 
cc-pair in the nuclear matter as exp{—{LA + Lb)/X)- 
Here La and Lb are the lengths for the trajectories of 
the cc in the projectile {A) and target (B) respectively. 
They are calculated from collision geometry by using 
the oft-used relations between mean transverse energy 
of the bin, Et, and the average impact parameter b. 

Figure [20] exhibits [ |57] the results of the NA50 col- 
laboration on J/^ cross section as a function of the 
transverse energy Ex in Pb-Pb collisions at ^/s ~ 17 
GeV. It is normalized to the Drell-Yan cross section 
in the mass range shown and _B^^ is the branching 
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Figure 21. J/V'-suppression in Pb-Pb colhsions at SPS 
as a function of the energy density e. Figure taken from 



fraction of Jf-ip in the dimuon channel. The full curve 
depicts the expected normal suppression as a function 
of Et, computed as explained above using the fitted 
J/V' cross section of 4.18 mb obtained from the NA50's 
own pA data. The dashed lines show the computed 
error bars on the expected suppression, and the inset 
shows the ratio of measured to the expected suppres- 
sion. Using the Bjorken formula in eq. ([4]), one obtains 
this ratio of the measured to the expected cross sec- 
tion ratio of the J/tp and the Drell-Yan as a function 
of the energy density in GeV/fm"^ units, as shown in 
Figure [2T1 taken from [f57]- One sees that the anoma- 
lous suppression, i.e., depletion of the measured cross 
section from that expected, sets in at an energy den- 
sity of about 2.5 GeV/fm'^, comparable to the expecta- 
tions from lattice QCD, as seen in Figure [J] A natural 
explanation of the anomalous suppression was, there- 
fore, the formation of quark-gluon plasma. Since the 
J/'(/'-production takes place both directly and through 
other charmonium states like Xc the slow fall-off with 
the energy density in Figure [^T] could be interpreted 
as gradual progress towards the full suppression. How- 
ever, one could also explain the anomalous suppression 
in alternative ways, using hadronic [ 58J or thermal [ 159] 
models. Since one expects the higher collision energy 
at RHIC to produce higher temperatures/energy den- 



18 



Bhalerao & Gavai 



sities, one expected a further stronger suppression at 
RHIC. Indeed, this seems to be true both in the quark- 
gluon plasma models as well as the alternatives, the 
difference between them being quantitative in nature. 
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Figure 23. Comparison of NA50 and PHENIX results 
on J/ f/j-suppression as a function of number of partici- 
pants. Figure taken from [I61j. 



Figure 22. J/-0-suppression in Au-Au collisions at 
PHENIX, BNL as a function of number of participants. 
Figure taken from [ [60] . 

The RHIC results [[SD], however, brought a big sur- 
prise by being different from any of those expectations. 
Analogous to the case of jet quenching in the previous 
section, the PHENIX collaboration at RHIC constructs 
the ratio Raa of the J/tp (differential) production cross 
section in AA collisions and the corresponding pp cross 
section weighted by the number of binary collisions. 
Figure [22] displays their results for Raa in Au-Au col- 
lisions at ^/s = 200 GeV. They show more suppression 
in the forward region {\y\ £ [1.2, 2.2], filled circles in the 
top panel), than the central {\y\ < 0.35, open circles in 
the top panel) for number of participants greater than 
100 (alternatively for large enough transverse energy 
Et). More importantly, a direct comparison [ '61] in 
Figure [23] clearly demonstrates that the PHENIX data 
in the central rapidity region are in very good agree- 
ment with the CERN NA50 resuhs ['37]. The trends for 
both the central region of the CERN and RHIC exper- 
iments, as seen in Figure [23] and the ratio of forward 
to the central rapidity region, as seen in the bottom 
panel of Figure [22l are against [ [61] the predictions of 
the models which successfully accounted for the NA50 



data. 

There have been some attempts to solve this J/ip- 
puzzle. As we saw in the Figure [TU] of section 2.4, the 
lattice QCD results suggest melting of the J/tp takes 
place at higher temperatures (> 2Tc) than predicted 
by simple models. A way to understand the results in 
Figure [23] could then suggest itself if the temperature 
reached at both the SPS and RHIC energy is ^ 2Tc. 
In that case, only Xc and ip' would have melted [ [62] , 
suppressing the corresponding decay J/ip's, and giving 
similar results for CERN and RHIC experiments. Since 
the temperature reached at LHC is expected to cross 
2Tc, a clear prediction of such a scenario would then 
be much more suppression for LHC than that in Figure 
[23l However, there are other scenarios, including ther- 
mal enhancement [ [63] arising due to recombination of 
the large number of thermal produced charm-anticharm 
quarks. These would predict an overall enhancement. 
In any case, J/-0-suppression could provide a lot of ex- 
citement again at LHC. 

3.5. Particle ratios & Bulk Properties 

A variety of hadrons are produced in an ultra- 
relativistic heavy-ion collision. They are identified and 
their relative yields measured; see Fig. [24] These 
hadron abundance ratios can be calculated in a sim- 
ple statistical model [ [Mj: It is assumed that these 
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particles emerge from a chemically equilibrated hadron 
gas characterized by a chemical potential {^i) for each 
hadron species and a common temperature (T). The 
number density rij of hadron of type i is then given 
by the standard Fermi-Dirac (+) or Bose-Einstein (— ) 
formulas 

^ , 1 

V (27r)3exp[(£;,-M.)/7^]±l' 

where di is the spin degeneracy. At chemical equi- 
librium, the chemical potential fii can be written as 

(3) (3) 

Mi — tJ'sBi — ^isSi — fiilj^ where Bi, Si and stand 
for the baryon number, the strangeness and the third 
component of the isospin quantum numbers, respec- 
tively, of the hadron of type i. The two unknown pa- 
rameters T and fiB are fitted to the data. This simple 
model has been quite successful in explaining the SPS 
and RHIC data; see Fig. [M] for SPS and a similar fig- 
ure in [[65] for RHIC. Note that even the multistrange 
particles seem to be consistent with the model. This 
suggests that they are produced in a partonic environ- 
ment rather than in a hadronic one. T = Tch is the 
chemical freezeout temperature. The fitted values are 



Tch = 170 MeV, HB 
Tch = 176 MeV, 
Tch = 177 MeV, ^ib 



= 270 MeV, (SPS), 

= 41 MeV, (RHIC 130 GeV), 

= 29 MeV, (RHIC 200 GeV). 



Note the trend of the chemical freezeout point to ap- 
proach the temperature axis of the QCD phase diagram 
as the collision energy is increased. Data obtained at 
the AGS and SIS energies are also consistent with this 
trend; see Fig. 1.3 in [[66]. For more recent fits to the 
statistical model, see [ [67] . 

4. Hydrodynamics 

Hydro plays a central role in modelling relativistic 
heavy-ion collisions: It is first used for the calculation 
of the pt spectra and the elliptic flow V2 ■ The resultant 
energy density or temperature profiles are then used in 
the calculations of jet quenching, J/ip melting, thermal 
photon and dilepton production, etc. 

Hydrodynamic framework consists of a set of coupled 
partial differential equations for energy density, number 
density, pressure, hydrodynamic four-velocity, etc. In 
addition, these equations also contain various transport 
coefficients and relaxation times. 

Hydro is a very powerful technique because given the 
initial conditions and the EOS it predicts the evolution 
of the matter. Its limitation is that it is applicable at 
or near (local) thermodynamic equilibrium only. 
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Figure 24. Comparison between the statistical model 
(horizontal bars) and experimental particle ratios 
(filled circles) measured at SPS CERN. From Braun- 
Munzinger et al. [ I64j . 



4.1. A Perfect Fluid? 

How robust is the claim of discovery of a perfect fluid 
at RHIC, or is there any need of the viscous hydro- 
dynamics for RHIC? A closer scrutiny shows that the 
claim is not really robust, and it is necessary to do vis- 
cous hydro calculations: 

• Agreement between data and ideal hydro is far 
from perfect. (Ideal) "hydro models seem to work for 
minimum-bias data but not for centrality-selected tt 
and p data" [ 68 . 

• Initial (and final) conditions for the hydrodynamic 
regime are uncertain. It is entirely possible that the 
ideal hydro mimics viscous hydro if the initial (and/or 
final) conditions are suitably tuned. Most ideal hydro 
calculations so far have been done with Glauber-type 
initial conditions. It has recently been realized that the 
CGC-type initial conditions yield higher eccentricity of 
the overlap zone [ i69j , and hence higher V2 ■ To push 
these results down to agree with data, some viscous cor- 
rections are needed. The same is true with fluctuations 
in the initial conditions [ TO] . Event-to-event fluctua- 
tions in nuclcon positions result in higher eccentricity 
and hence higher V2 [ 171] . 

• Some V2 may build up during the pre-equilibrium (i.e., 
pre-hydro) regime. Success of ideal hydro may be due 
to the neglect of this contribution to V2 in most calcu- 
lations [[72]. 

• For realistic light quark masses, the deconfinement 
transition is known to be a smooth crossover. However, 
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it seems that the ideal hydro calculations need a first- 
order transition for a best fit to the data [ 73 • 

• The shear viscosity to entropy density ratio {rj/s) may 
be small in the transition region. But there are indi- 
cations that the bulk viscosity to entropy density ratio 
(C/s) may be rising dramatically near Tc [IZ^- If this 
result holds, QGP discovered at RHIC cannot be called 
a perfect fluid. 

• It is known that for helium, water, nitrogen, rj/s at 
constant pressure plotted as a function of temperature, 
exhibits a minimum with a cusp-like behaviour at the 
critical point; see Fig. [23 There are indications that 
the QCD matter too shows similar trends. Viscous hy- 
dro calculations of the QCD matter would allow us to 
extract rj/s from data and might help us pinpoint the 
location of the QCD critical point [I75j. 
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Figure 25. Each curve is at a fixed pressure. Solid: 
below the critical pressure Pc, dotted: at Pc, dashed: 
above Pc- From [[75]. 

• If the inequality rj/s > l/An obtained [ 16J from the 
AdS/CFT duality is applicable to QCD, then also vis- 
cous hydro calculations become necessary. 

• Assume a quasiparticle picture. Quantum mechanical 
uncertainty principle tells us that the mean-free path 
(A) cannot be less than the inverse of the typical mo- 
mentum of the quanta. It also makes no sense to have 
a mean-free path smaller than the interparticle spacing 
[ 176] . Since rj on X, cannot vanish. 

• Finally, to claim success for ideal hydro, one should 
calculate viscous corrections and show explicitly that 
they are indeed small. 



4.2. Relativistic Dissipative Hydro — a Brief 
History 

Relativistic version of the Navier-Stokes equation was 
obtained by Eckart [ 77J, and by Landau and Lifshitz 
[ [78] . This is called the standard or the first-order for- 
malism because terms only up to first order in dissipa- 
tive quantities are retained in the entropy four-current. 
(The Euler's equation constitutes the zcroth-order for- 
malism.) However, it was soon realized that this for- 
malism suffers from the following problems: 

• Acausality: Equations are parabolic and they result 
in super- luminal propagation of signals [ 79, 80J. 

• Instability: Equilibrium states are unstable under 
small perturbations for a moving fluid ['8T\. This makes 
it difficult to perform controlled numerical simulations. 

• Lack of relativistic covariance: This problem is re- 
lated to the previous one. First-order theories look co- 
variant, but they are not. 

A causal dissipative formalism was developed by 
Miiller [ |79] , and Israel and Stewart [ [80] , in the non- 
relativistic and relativistic sectors, respectively. It is 
also called a second-order formalism because the en- 
tropy four-current now contains terms up to second or- 
der in dissipative quantities. The resulting hydrody- 
namic equations are hyperbolic. Application of causal 
dissipative hydro to relativistic heavy-ion collisions was 
pioneered by Muronga [ 82 . Since then many others 
have contributed to this effort. We shall describe some 
of them in subsection 4.4. 

Recent years have witnessed intense activity in the 
area of causal hydro of gauge theory plasmas from 
AdS/CFT duahty; for reviews see [[83|. 

4.3. Basic Idea of Causal Dissipative Hydro 

Before we discuss hydrodynamics, let us first con- 
sider a simpler example of diffusion. Consider a fluid 
in equilibrium with a uniform density p. If the fiuid is 
perturbed such that the density is no longer uniform, it 
responds by setting up currents which tend to restore 
the equilibrium. In the linear response theory, the in- 
duced current Ji is simply proportional to the gradient 
of p (Fick's law.): 

J, = -Dd^p, (8) 

where D is the diffusion coefficient. D is an example of 
a transport coefficient. Transport coefficients play an 
important role in the study of relaxation phenomena in 
non-equilibrium statistical mechanics or fiuid dynam- 
ics. Equation ([S]) connects the applied force {~dip) 
with the flux (Ji). Such equations are called constitu- 
tive equations because they describe a physical prop- 
erty of the material. (The familiar Ohm's law J = ctE 



Heavy Ions at LHC: A Quest for Quark-Gluon Plasma 



21 



is another example of this.) In addition to eq. ([5]), we 
also have the usual current conservation equation 

9p = 0. (9) 

If D is constant, elimination of Ji gives 

dop - Ddfp = 0. 

This is the diffusion equation. It is parabolic. Its solu- 
tion is 

p ~ exp(— a; 

It is easy to see that the solution violates causality: Ini- 
tially (i.e., in the limit t — > 0), this is the Dirac delta 
function. But at any finite time, howsoever small, it is 
nonzero everywhere, even outside the lightcone. Now 
eq. ([9]) cannot be wrong. So to restore causality the 
constitutive equation ([5]) which anyway was a hypoth- 
esis, is replaced by 

TjdoJ^ + J, = -Dd,p, (10) 

where tj is a parameter with dimensions of time. In 
eq. ([U, if the force vanishes, the flux vanishes instanta- 
neously without any time lag. In contrast, in eq. (jlOp 
the flux relaxes to zero exponentially, tj is called the 
relaxation time. The new diffusion equation is 

Tjdlp + d^p~Ddjp^Q. 

This equation is hyperbolic and is called the Tele- 
graphist's equation [ 84 . If = D/tj < 1, causality 
is restored. 

Now consider hydrodynamics. The conservation and 
constitutive equations are 

d^T''-' - 0, 

2 

Tij = P5ij - ri{diUj + djUi- -SijdkUk) 
- C^ijdkUk. 

Here T^'^ is the energy-momentum or stress-energy ten- 
sor, P is the equilibrium pressure, and 77 and ^ are 
the coefficients of shear and bulk viscosity, respectively. 
Tensor decomposition is now more complicated. But 
the basic idea remains the same. Causality is restored 
by introducing higher-order terms in the gradient ex- 
pansion. This forces introduction of a new set of trans- 
port coefhcients, e.g., r^r and rn which are relaxation 
times corresponding to shear and bulk viscosities. They 
are important at early times or for a rapidly evolving 
fluid. For details, see e.g. [|82j. 



4.4. Recent Results from Causal Viscous Hydro 

The Israel-Stewart formulation [ 80J of the causal dis- 
sipative hydro is commonly used for numerical applica- 
tions. However, it is not the only causal formulation 
available. There are others such as Miiller's theory [ 
[75] . Carter's theory [[SS], Ottinger-Grmela formulation 
[ [55] , memory function method of Koide et al. [ [57] , 
etc. 

We have already mentioned the early work by 
Muronga [[82]. Since then several authors have studied 
various aspects of the causal viscous hydro. We now 
describe briefly only a few of the most recent of these 
papers. This will also give the reader a feel for the 
complexities of these calculations and the uncertainties 
therein. (Other very recent papers which we shall not 
describe are listed in [ 8.8].) 

Romatschke and Romatschke [ ^531 used the Israel- 
Stewart theory. They assumed longitudinal boost in- 
variance and used Glauber- type initial conditions. The 
initial shear pressure tensor tt'^" was assumed to be 
zero, ri/s was treated as a fixed number independent 
of temperature. The bulk viscosity was ignored. For 
the EOS they used the semirealistic result of Laine and 
Schroder [ 90J, and calculated the elliptic flow V2- Their 
conclusion was that pT-integrated V2 is consistent with 
rj/s up to 0.16; see Fig. 1261 However, the minimum-bias 
V2{pt) favoured rj/s < l/Air violating the KSS bound [ 
[IS]; see Fig. [571 




Figure 26. Au-Au, 200 GeV, pj-integrated V2 for 
charged particles vs number of participant nucle- 
ons. PHOBOS: 90% confidence level systematic errors. 
From [Mi- 
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whereas C/s decreases with T. This is easy to under- 
stand because QCD becomes conformally symmetric at 
high temperatures. 

In the deconfinement transition region the conformal 
symmetry is badly broken, and there is no reason to 
expect the bulk viscosity to be negligible. Extracting 
( for temperatures in this region from lattice QCD is 
difficult; see section 2.4. However, some preliminary 
results are now available, and they indicate a dramatic 
rise of C/s asT^Tc [|74]. 



Figure 27. Au-Au, 200 GeV, minimum-bias V2{pt) for 
charged particles. STAR: only statistical errors. From 



Dusling and Teaney [W^. used the Ottinger-Grmela 
formalism of causal viscous hydro. They assumed longi- 
tudinal boost invariance and used Glauber-type initial 
conditions. The initial shear pressure tensor tt'^ was 
taken to be rj(d'''u^) as in the Navier-Stokes theory. 
77/s was treated as a fixed number independent of tem- 
perature. The bulk viscosity was ignored. The EOS 
used by them was simply p — e/i without any phase 
transition. Their conclusion was that if the effects of 
viscosity are included in the evolution equations but 
not in the freezeout, then the V2 is affected only mod- 
estly. If, however, they are included at both the places, 
then V2 is significantly reduced at large pt- 

Why does the shear viscosity suppress 1^2 (pt)? Shear 
viscosity represents a frictional force proportional to 
velocity. For an in-plane elliptic flow, the in-plane flow 
velocity is higher than that out of plane. So the in- 
plane frictional force is stronger. This tends to reduce 
the flow anisotropy and hence V2{pt)- 

Calculations described above include the shear vis- 
cosity in some approximation, but ignore the bulk vis- 
cosity completely. What do we know about the bulk 
viscosity of the strongly interacting matter? In the 
high-temperature limit, pQCD calculations [ [22] give 
the following results for the shear and bulk viscosity 
coefficients 



rp3 



and C, ' 



In Us 



As T increases, both rj and C, increase. However, the 
ratio C/rj decreases showing the reduced importance of 
the bulk viscosity at high T. Also note that the en- 
tropy density s ^ , and hence rj/s increases with T, 




Figure 28. Bulk viscosity based on lattice data, loq = 
0.5,1,1.5 GeV (top to bottom) is the scale at which 
pQCD is applicable. From [[74]. 



Taking these results at their face value, Fries et al. [ 
[Sni have studied the effect of inclusion of the bulk vis- 
cosity in the hydro equations. They studied ID expan- 
sion of the fluid assuming longitudinal boost invariance. 
-q/s was held fixed at 1/Att. A realistic EOS based on 
the lattice results of Cheng et al. [[12] was used. Vari- 
ous initial conditions were tried. They concluded that 
(a) Large bulk viscosities around Tc lead to sizeable 
deviations from equilibrium throughout the entire life- 
time of QGP. (b) Bulk viscosities just slightly larger 
than currently favoured could easily lead to breakdown 
of hydro around Tc- (c) The decreased pressure should 
slow down the expansion and increase the time spent by 
the fiuid in the vicinity of the phase transition, (d) The 
amount of entropy produced through bulk stress around 
Tc is smaller than that produced by shear stress at ear- 
lier times. Hence no large increase of the final particle 
multiplicity is expected. 
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4.5. What Remains to be Done? 

• Bulk as well as shear viscosity (together with temper- 
ature dependence of C/s and rj/s) needs to be incorpo- 
rated. 

• Can causal viscous hydro with CGC-type initial con- 
ditions reproduce dN/dy, (pr) and V2 data? If so, what 
are the extracted C/s, ?7/s? 

• Causal viscous hydro + hadronic cascade is not done 
yet. 

• There are issues related to the hydro formalism itself. 
For example, Baier et al. [[94] have recently shown that 
the Miiller and Israel-Stewart theories do not contain 
all allowed second-order terms. 

• Present uncertainties in the hydro calculations limit 
the accuracy with which conclusions can be drawn. A 
coherent, sustained collaboration of experts in all stages 
of heavy- ion collisions is needed for a detailed, quantita- 
tive analysis of experimental data and theoretical mod- 
els. Various numerical codes need to be compared with 
each other. To that end a new Theory- Experiment Col- 
laboration for Hot QCD Matter (TECHQM) has been 
initiated. For details, see [I95j. 

5. Predictions for LHC 

Pb-Pb collisions at ^snn — 5.5 TeV is an important 
part of the LHC experimental program. 5.5 TeV rep- 
resents about 30-fold increase in the CM energy com- 
pared to the maximum energy explored at RHIC which 
in turn was about 10 times higher than that at SPS. 
Measurements on pp collisions as well as collisions of 
p, d, light ions with Pb will provide important bench- 
marks. 

Among the experiments at LHC, CMS and ATLAS 
are primarily particle physics experiments/detectors, 
but they will study the physics of heavy-ion collisions 
too. ALICE (A Large Ion Collider Experiment), on the 
other hand, is a dedicated heavy-ion collision exper- 
iment. Physicists from several Indian universities and 
institutions have contributed in a big way to the ALICE 
collaboration. They are responsible for, among other 
things, the designing, testing, installation and mainte- 
nance of the Photon Multiplicity Detector (PMD) in 
ALICE and future upgrades of it. PMD is a preshower 
detector with fine granularity, full azimuthal coverage 
and one unit of pseudo-rapidity coverage. It will be 
used to measure the multiplicity, spatial distribution 
and correlations of produced photons on an event-by- 
event basis. Since photons escape the quark-gluon 
plasma without interactions, such measurements can 
potentially provide a cleaner glimpse of the early QGP 
phase. The Indian community has also made signifi- 



can contributions to the muon spectrometer of ALICE. 
The spectromemter will be useful in the investigations 
of the J/^jj and other quarkonia, discussed in subsection 
3.4. These particles are detected via their dimuon decay 
channel. The muon tracks will be found with an accu- 
racy of better than one-tenth of a millimeter, thanks 
to the state-of-the-art readout electronics, known as 
MANAS, which was developed indigenously. ALICE 
has decided to use a Grid environment for their com- 
puting needs. India is a signatory to the Worldwide 
LHC Computing Grid and some of the Department of 
Atomic Energy installations are designated as Tier-II 
centers for this purpose. 

A workshop was organized in 2007 at CERN in order 
to collect all the existing predictions for heavy-ion col- 
lisions at LHC. The proceedings [[M] provide a broad 
overview of the field. Here we shall only present a few 
glimpses of what may be in store at LHC. 
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Figure 29. Charged-particle rapidity density per partic- 
ipant pair as a function of center-of-mass energy for AA 
and pp collisions. Dashed line: a fit linear in ln(-y/s). 
Dotted curve: a fit quadratic in In(Y^), Long-dashed 
curve: based on the saturation model of [[97]. From [ 
[66]. 

One of the first and easiest measurements at ALICE 
would be that of the charged-particle multiplicity in the 
mid-rapidity region. Particle production models and 
simple fits which are in agreement with the AGS, SPS, 
and RHIC data on this quantity differ substantially 
from each other when extrapolated to the LHC energy, 
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Figure 30. Pseudorapidity-azimuthal angle plot of Pb- 
Pb event at LHC energy with two 100 GeV jets gen- 
erated with HIJING and PYTHIA event generators. 
From [[98]. 



as shown in Fig. [291 Thus this simple "first-day" mea- 
surement will test our understanding of the physics of 
multiparticle production. The charged-particle multi- 
plicity provides a handle on the initial entropy produc- 
tion; the latter quantity is a necessary input in the hy- 
drodynamic evolution of the produced matter. 

Another relatively simple measurement at ALICE 
would be that of the elliptic flow V2 which has played 
a crucial role at RHIC (sec. 3.2). The initial energy 
density (eq. (|4|)) as well as the QGP lifetime arc pre- 
dicted to be higher at LHC than those at RHIC. This 
is expected to raise the value of V2{pt)- On the other 
hand, the increased radial flow at LHC is expected to 
lower it. (Recall the discussion on mass ordering in sec. 
3.2.) The net effect on V2{pt) depends on the mass of 
the hadron: Minimum-bias V2{pt) for pions (protons) 
is expected to be higher (lower) at LHC than at RHIC, 
at low pt\ see Eskola et al. in [[96j. Prediction by 
Kestin and Heinz is that V2(j>t) at a fixed impact pa- 
rameter will be smaller at LHC than at RHIC, for pions 
as well as protons [[96]. However, pT-integrated elliptic 
flow is expected to be higher for all hadrons due to the 
increased relative weight at large values of px ■ 

In sec. 3.5 we have quoted the values of Teh and fiB 
for the SPS and RHIC energies. The latest predictions 
for LHC are T^h = 161±4 MeV and ^lB = O.StJ g MeV 

m- 

Hard processes: Cross sections for the production of 
heavy flavours, CTcc and fTf,6' are expected to be about 
10 and 100 times larger at LHC than at RHIC. Cross 



sections for the production of jets with transverse en- 
ergy in excess of 100 GeV are expected to be several 
orders of magnitude higher. Jet-photon events will also 
be abundant. Figure [30] displays the capability of AL- 
ICE to reconstruct the high-energy jets at LHC in spite 
of the large soft-hadron background. Thus it would be 
possible to make detailed differential studies of heavy- 
quarkonium production, open-charm and open-beauty 
production, jet quenching, etc. at LHC [[9^. It will 
also be possible to study quark mass dependence and 
colour charge dependence of the energy loss of a parton 
as it traverses the medium. 

Thus LHC promises to be a valuable tool to test 
our models of ultrarelativistic heavy-ion collisions and 
deepen our understanding of QCD. For details, see [ 
[99]. 
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